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We consider a simplified model for self-induced flavor conversions of a dense neutrino gas in two 
dimensions, showing new solutions that spontaneously break the spatial symmetries of the initial 
conditions. As a result of the symmetry breaking induced by the neutrino-neutrino interactions, 
the coherent behavior of the neutrino gas becomes unstable. This instability produces large spatial 
variations in the flavor content of the ensemble. Furthermore, it also leads to the creation of domains 
of different net lepton number flux. The transition of the neutrino gas from a coherent to incoherent 
behavior shows an intriguing analogy with a streaming flow changing from laminar to turbulent 
regime. These finding would be relevant for the self-induced conversions of neutrinos streaming-off 
a supernova core. 

PACS numbers: 14.60.Pq, 97.60.Bw 


I. INTRODUCTION 

Neutrino-neutrino interactions in dense neutrino gases 
provide a refractive term leading to a non-linear feed¬ 
back in the flavor evolution, when the neutrino interac¬ 
tion potential /i ^ \/2GFnu is comparable or larger than 
the neutrino vacuum oscillation frequency w = /2E. 

This effect can lead to a collective behavior of the dense 
neutrino ensemble (see [I| for a recent review). Remark¬ 
ably, almost a decade ago it was realized that neutrino- 
neutrino interactions dominate the flavor evolution in the 
deepest core-collapse supernova ^N) regions producing 
self-induced flavor conversions [2 t 3| • The most important 
observable consequence of this type of flavor transitions 
would be the swap of the SN Vg and Ue spectra with that 
of non-electron flavor and Ux in certain energy inter¬ 
vals 

Characterizing the SN neutrino flavor dynamics 
amounts to follow the spatial evolution of the neutrino 
fluxes. For a stationary neutrino emission, the kinetic 
equations of the u space-dependent occupation numbers 
Pp.x with momentum p at position x are 00 

Vp ' Vx Pp,X — ^[f^p.X; f?p,x] 7 (1) 

with the Liouville operator in the left-hand side. Ne¬ 
glecting external forces and an explicit time dependence 
of the occupation numbers, this operator represents the 
drift term proportional to the neutrino velocity Vp, due 
to particle free streaming. On the right-hand-side the 
matrix rip^x is the full Hamiltonian containing the vac¬ 
uum, matter and self-interaction terms. In particular, in 
non-isotropic neutrino gases, like the case of neutrinos 
streaming-off a SN core, the neutrino-neutrino interac¬ 
tion term contains multi-angle effects since the current- 
current nature of the low-energy weak interactions intro¬ 


duces an angle dependent term (1 — Vp • Vq) between two 
interacting neutrino modes 0 0 . This term produces a 
net current so that test neutrinos moving in different di¬ 
rections would experience different refractive index. This 
in some cases challenges the collective behavior of the fla¬ 
vor evolution observed in an isotropic case, leading to fla¬ 
vor decoherence im. Multi-angle effects can also lead 
to a trajectory-dependent matter term, which if strong 
enough suppresses the self-induced conversions 

The multi-angle flavor evolution described by the par¬ 
tial differential equations dll) has never been solved 
till now in its full complexity. Instead, numerical ap¬ 
proaches have been typically based on the so-called “bulb 
model” iiini, where it is assumed a spherical sym¬ 
metry about the center of the SN and azimuthal sym¬ 
metry about any radial direction. In this limit Eq. 0 
reduces to an ordinary differential equation problem, 
projecting the evolution along the radial direction, i.e. 
Vp • Vx — Vr d/dr. 

Attempts to go beyond the bulb model have been pro¬ 
posed. For example, in 0 it was shown that assuming 
that the neutrino ensemble displays self-maintained co¬ 
herence, the problem for generic geometries can be re¬ 
duced to a one-dimensional case along the streamlines 
of the overall neutrino flux. However, the existence of a 
self-consistent coherent solution does not imply its sta¬ 
bility. Indeed, it has been recently realized that insta¬ 
bilities may grow once one relaxes some symmetries of 
the bulb model, since neutrino-neutrino interactions can 
lead to spontaneous symmetry breaking effects. In 0 it 
was shown that perturbing the symmetry of the initial 
neutrino angular distributions one would find a speed¬ 
up in the flavor exchange. Moreover, removing the as¬ 
sumption of axial symmetry in the v multi-angle term, 
a multi-azimuthal-angle instability has been discovered. 
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even assuming a perfect spherically symmetric v emis¬ 
sion [iM 3 - Furthermore, also space and time homo¬ 
geneity can be broken in a dense neutrino gas, so that it 
is not guaranteed that a quasi-stationary neutrino emis¬ 
sion would lead to a stationary solution Finally, 

with a simple toy model it has been recently shown, by 
means of a stability analysis of the linearized equations of 
motion, that self-induced oscillations can spontaneously 
break spatial symmetries [2l[. All these findings suggest 
that the validity of the bulb model should be critically 
reconsidered and that a self-consistent solution of the SN 
neutrino flavor evolution can only be achieved by solving 
the complete multi-dimensional problem of Eq. (HI)- 
As a further step in clarifying this issue we consider 
here the two-dimensional toy model discussed in [^ . 
i.e. monochromatic neutrino streams emitted in a sta¬ 
tionary way in two directions from an infinite boundary 
plane at z = 0 with periodic conditions on x and trans¬ 
lation invariance along the y direction. We assume a 
small perturbation for the initial symmetries of the fla¬ 
vor content in both the two emission modes and along 
the boundary in the x direction. Despite of its sim¬ 
plicity, this model with perturbed symmetries exhibits 
a rich phenomenology. Indeed, solving numerically the 
flavor evolution we find that the initial small perturba¬ 
tions are amplified by neutrino interactions, leading to 
non-trivial two-dimensional structures in the oscillation 
pattern that exhibits large space fluctuations. Any par¬ 
ticular solution which is required to satisfy initial sym¬ 
metries at the boundary is thus, unstable, confirming the 
results of [2l[. We will show that the neutrino flux tran¬ 
sition from a coherent to incoherent behavior resembles 
quite closely the transition from a laminar to turbulent 
regime of a streaming flow. This issue will be studied in 
more details elsewhere. Finally, also the lepton number 
flux would develop a domain structure. 


II. A TWO-DIMENSIONAL MODEL 

As in , neutrinos are emitted from an infinite plane 
at z = 0 , in two directions 

W = (C = L,R) , (2) 

where 0 < 1)2 < 1 and vr = —vl = 1/I — The 
flavor evolution occurs for z > 0 in the two-dimensional 
plane spanned by the x and z coordinates. We consider 
a two-flavor {ve,Vx) neutrino ensemble. Expanding all 
quantities of Eq. (HD in terms of Pauli matrices ct, one 
gets the equations for the L mode 

'Vl-^J^l{x,z) = [-\-ujB + fiDR{x,z)] X Pl{x,z) , 

Vl • VxPila;, z) = [-wB-f/iDfl(a;, z)] x Pl(x, z) ,(3) 

with u! = Am'^/2E the vacuum oscillation frequency. 
The equations from the R mode can be obtained from 
the previous ones using the L R symmetry. 


The sans-serif in Eq. indicates three-dimensional 
vectors in flavor space, e.g. B = , B^). The dif¬ 

ferential operator at left-hand-side reads 

Vl ■ Vx = vlOx -b Vzdz . (4) 

The Pl,r (Pl,/?) functions are the neutrino (antineu¬ 
trino) polarization vectors in flavor space for the L, R 
modes, defined from the occupation numbers as in [ll| . 
i-e. Ql,r = -b -_Fy^)PL.,R ■ cr] and 

Ql,r = l/ 2 [^'pe + Ppx + {Fi>e - - Fpx ) PL,fl • cr], where Fp 
are the total neutrino fluxes of the different species. We 
define as usual B^r = Pl,r — Pr.r- The unit vector B 
points in the mass eigenstate direction in flavor space, 
such that B • 63 = — cos0, where 9 is the vacuum mixing 
angle. For simplicity we neglect a possible matter effect, 
assuming that its only role would be to reduce the ef¬ 
fective in-medium mixing angle, 6 1 [3. Finally, the 

neutrino self-interaction term is given by [ll| 

y = y/2GF[Fl - f 0J(1 - • v^) , (5) 

expressed in terms of the antineutrino fluxes Fp at 
the boundary. Since in our model the neutrino tra¬ 
jectories intersect at a fixed angle, the previous equa¬ 
tion would correspond to the “single-angle” approxima¬ 
tion [^. Therefore, multi-angle effects do not play a role 
in this case. 

One can define a conserved “lepton current” = 
(Lo,L) whose components are (see [^ 1 

Lq = Dl • B + D/j • B , (6) 

L = Vi(Di.B)+Vfl(Dfl.B) , (7) 

where L is a two-dimensional vector in the {x, z) plane, 
o —3 

and Dl,r • B ~ F£ ^ — F^ In the following we will 
normalize to the number of modes Nr^r = 2. It is 
straightforward to show from Eq. that L satisfies a 
continuity equation 

^iLo + Vx • L = Vx • L = 0 , (8) 

where first equality follows from the assumption of a sta¬ 
tionary solution. Eq. ([HD generalizes the lepton-number 
conservation law of the one dimensional case Q. 

As first exploited in (see also mi) in the context of 
multi-dimensional neutrino oscillations, the partial differ¬ 
ential equation problem like of Eq. ([31 can be reduced to 
a tower of ordinary differential equations for the Eourier 
modes defined as 

/ -I-00 

dx PL(R)ix,z)e~'^^^ , (9) 

-00 

and similarly for the antineutrino polarization vectors 
Pl,r- In the following, for simplicity we will consider 
a monochromatic perturbation in the x direction for 
the neutrino polarization vectors at the boundary at 
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FIG. 1: Three dimensional evolution of the flavor content 
flee in the x-z plane, and its map on the bottom plane. 


FIG. 2: Map of the lepton number Lq in the x-z plane. 


z = 0. Since we start with a pure flavor state, we as¬ 
sume Pi f^{x, 0) = Pi j^{x, 0) = 0 and 

0 ) = 0)) + ecos(fcox) , (10) 

where this latter component of the polarization vector is 
proportional to the flavor content of the ensemble. The 
function {P^ nix.O)) indicates the unpertubed value of 
the polarization vectors, while fco is the wave-number of 
the perturbation and e /r, u; its small amplitude. It is 
easy to see that in this case, only higher harmonics of the 
fundamental mode with = nk^ are excited. Defining 
Pl,™ = fcoPL.fc„/(27r), from Eq. ([3]) one obtains 


j PL,n('^) — L,n “b ^ 

dz 


+ 00 

j=-oo 


An analogous set of coupled ordinary differential equa¬ 
tions can be written for the R mode and for the antineu¬ 
trino polarization vectors. It is enough to follow the evo¬ 
lution for positive modes, n > 0, since the PL,Rix, z) and 
Pl.r{x, z) are real functions and therefore 


P(L,R),n = PiL,R),-n ■ (12) 

Once the evolution of the harmonic modes is obtained 
from Eq. (113, the polarization vector in configuration 
space can be obtained by inverse Fourier transform. 


III. NUMERICAL EXAMPLES 

To illustrate the behavior of the self-induced flavor con¬ 
versions in our two-dimensional toy model, we consider a 
gas initially composed by only z/g and Pe- We normalize 
the polarization vectors to the z/g number density and we 
assume an excess of z^g over i/g, i.e. F^/Pp^ = I -|- a. 
If the translational symmetry is assumed, i.e. e = 0 in 


Eq. (flOl) . and the L and R modes are prepared identically, 
it is well known that the system is stable in normal mass 
hierarchy {Am? > 0), while in inverted mass hierarchy 
{Am? < 0) it exhibits a bimodal instability and behaves 
as a flavor pendulum, leading to periodic pair conversions 
zzgZ/g ^ VxPxjL that conserve the lepton number Lq = a 
of Eq. ITSl) Q . More recently, it has been shown that if 
the L O i? symmetry is perturbed the system becomes 
unstable also in normal hierarchy, with a similar pendu¬ 
lum behavior [I^ . Our further step is to perturb also the 
translational symmetry at the boundary. 

We fix in Eq. dS]) ^ = 10, w = I, 0 = 10“^, and 
we consider the normal mass hierarchy case (the result 
would be similar for the inverted mass hierarchy). We 
take as asymmetry parameter a = 0.3. We assume 
vr = Vz = ■\/2/2. Perturbation in the L R symmetry 
are introduced by a 1% difference in the initial conditions 
between these modes. Furthermore, to perturb the trans¬ 
lational symmetry along the x direction, we assume in 
Eq. (fTOjl e = 0.01, and we take as perturbation frequency 
fco = where the square-root expression is the 

proper frequency of the unperturbed flavor pendulum Q . 
Correspondingly, ^(0) = Pr i(0) = e (and analogously 
for antineutrinos), while the higher order harmonics are 
initially vanishing. We follow the evolution of the first 
N = 600 Fourier modes. 

In Fig. 1 we show the flavor evolution of Ue flavor con¬ 
tent Qee{x,z), in the x-z plane summed over the L and 
R contributions and normalized to = 2. We see 

that till z ~ 2.5 the flavor evolution shows a transla¬ 
tional symmetry, being uniform in x direction. The fla¬ 
vor content presents the known pendular nutations in z 
direction, observed in the one dimensional evolution [i3| . 
Contrarily to what previously assumed, this coherent be¬ 
havior is not stable. Indeed, at z ~ 2.5 the transla¬ 
tional symmetry is broken. Then, g^e develops large 
variations in the x direction at length scales that be¬ 
comes smaller and smaller increasing the distance from 
the boundary. Furthermore, also the L R symme¬ 
try is broken, since Pl{x,z) and Pl{x,z) present large 
differences (not shown), in agreement with the stability 
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analysis of [2l|. The two-dimensional map of the func¬ 
tion Qee{x, z) further clarifies the effect of the breaking of 
the translational symmetry. Till the symmetry is unbro¬ 
ken all the neutrinos oscillate in phase and the surfaces of 
equal phase are planes parallel to the radiating surface at 
z = 0. Then, when the instabilities develop these planes 
of common phase are broken and the coherent behavior 
of the oscillations is lost. 

In Fig. 2 we represent the lepton density Lq of Eq. ca) 
in the x-z plane. Notice that as soon as the translational 
instability develops, lepton number shows a non trivial 
domain structure and that self-induced conversions lead 
to large space variations of the initial asymmetry a. 

In order to understand the origin of this flavor dy¬ 
namics, in Fig. 3 we show a contour plot representing 
the evolution of the different Fourier modes |P_r,ti(-z)| 
(in logarithmic scale) in the plane of n-z. The behav¬ 
ior of \PL,n{z)\ would be similar (not shown). We realize 
that the breaking of the translational symmetry corre¬ 
sponds to the growth of the n > 0 modes occurring at 
z > 2. This dynamics can be seen as a cascade pro¬ 
cess in the Fourier space, where a flavor wave caused 
by the flavor pendulum diffuses to higher harmonics (i.e. 
to smaller scales) as soon as the Fourier modes are ex¬ 
cited by the non-linear interactions between the different 
modes. Note the analogy of this process with the multi¬ 
angle decoherence associated with a diffusion of excita¬ 
tions in the multipole space Correspondingly, in the 
flavor evolution one observes the developments of spa¬ 
tial variations in the x directions at smaller and smaller 
scales. In the example we are studying, at z = 12 about 
the hrst 300 harmonics are significantly excited. Indeed, 
the number of harmonics that one follows determine the 
range of validity of the numerical simulation. We checked 
that the number of excited harmonics is sensitive to the 
neutrino-neutrino interaction potential since this fac¬ 
tor determines the strength of the terms responsible for 
the growth of the modes in the second line of Eq. (HZl). 
The growth of the harmonics is also enhanced with the 
initial flavor asymmetry a. Indeed, in the sum at right- 
hand-side of Eq. (El) the term increases with 

the initial flavor asymmetry, pumping the higher order 
harmonics. We comment that in realistic cases (e.g. for 
supernova neutrinos) one typically has a declining neu¬ 
trino potential /i. In this situation one would expect a 
cut-off in the number of excited modes in function of the 
distance from the boundary. 

The multi-angle extension of the results discussed here 
are presented in Appendix A. 

The behavior of the neutrino gas in our model has a 
nice analogy with the transition between laminar and 
turbulent behavior of a streaming fluid (see, e.g., [^b 
In this respect, it is useful to define an average neutrino 
velocity for the neutrino fluid at a point x = (x, z). Con¬ 
sidering for example the z/g flux, one has 


(Ve)x 


Qee,L^L + Qee,R'^ R 
Qee,L Qee,R 


(13) 



n 


FIG. 3: Contour plot of the first 600 Fourier mode |Pfl,n| (in 
logarithmic scale) in the plane n-z. 



FIG. 4: Streamlines of the flux (in vertical direction). 


Till the L ^ R symmetry is unbroken,we have (ve)x — 
Vz- Then, when a L-R asymmetry with variations in the 
X direction is produced, the average velocity starts to 
acquire a transverse component in the x direction. The 
“streamlines” of the neutrino flux are the solutions of 


^ _ (Ve)x 
ds |(Ve)x| 


(14) 


where s is a parameter along the line (see H). 

In Fig. 4 we show the streamlines defined in Eq. CK in 
vertical direction). One clearly sees that the transition 
between the coherent to incoherent flavor behavior ob¬ 
served in Eig. 1 corresponds to the change from a laminar 
to a turbulent regime. As soon as the translational sym¬ 
metry is broken, the streamlines become irregular and 
are no longer parallel to the z directions. Moreover, they 
exhibit large variations, as in the turbulent motion in a 
fluid, and tend to converge in preferred directions. 

The behavior we found has a fascinating similarity with 
the non-linear instabilities of fluid flows, described by the 
Navier-Stokes equation, see e.g. [3,13, as shown in the 
Appendix B. 
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IV. CONCLUSIONS 

We presented the results of a study of self-induced fla¬ 
vor conversions in a simple two-dimensional model. As 
predicted from the stability analysis performed in 
we find that the self-interacting neutrino gas can break 
the spatial symmetries of the initial conditions and the 
neutrino flavor content achieves large space variations. 
This implies that the coherent behavior of the neutrino 
gas is unstable under small spatial inhomogeneities. This 
process also leads to the formation of domains of different 
net lepton number Lq. 

There is a close analogy between the symmetries of 
the planar two-dimensional model considered in our work 
and the spherical SN case. Indeed, the translational sym¬ 
metry in our model corresponds to the spherical symme¬ 
try of the bulb-model and the L-R symmetry is equiva¬ 
lent to the axial symmetry in the spherical case. Never¬ 
theless, our toy model is much simpler than any realistic 
SN neutrino case, since we are neglecting multi-angle ef¬ 
fects, continuous energy spectra, ordinary matter effects 
and declining neutrino densities. All these effects would 
add additional complications and numerical challenges. 
However, if our results would apply also to the SN case 
this would radically change the current description of the 
self-induced flavor conversions, and would have interest¬ 
ing phenomenological consequences, like e.g., the gener¬ 
ation of a self-induced direction-dependent asymmetry 
in the lepton number flux. These intriguing possibili¬ 
ties call for further efforts to go beyond the bulb model. 
This challenge would be quite demanding in terms of new 
computig time and/or novel approaches. 
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Equations of motion 
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Va ■ VxPa(a;, 2;) = [-fwB 
N 

-K - Vg ■ V/3)D/3] X Pa{x,z), 

13=1 


Va ■ VxPa(a;, 2:) = [-wB 
N 

-t- - Vg • V/3)D;g] X Pa(a;, z) . 

f3=l 


The lepton number in this case is defined as 


A.M. acknowledges useful discussions with Eligio Lisi, 
Antonio Marrone and Georg Raffelt. A.M. also thanks 
Lun-Shin Yao for useful discussions on non-linear insta¬ 
bilities in fluids. The work of A.M. was supported by the 
German Science Foundation (DFG) within the Collab¬ 
orative Research Center 676 “Particles, Strings and the 
Early Universe” in Germany, and by the Italian Ministero 
deiristruzione, Universita e Ricerca (MIUR) and Istituto 
Nazionale di Fisica Nucleare (INFN) through the “The¬ 
oretical Astroparticle Physics” projects in Italy. G.M. is 
supported by INFN I.S. TASP. N.S. acknowledges sup¬ 
port from the European Union FP7 ITN INVISIBLES 
(Marie Curie Actions, PITN- GA-20II- 289442). 


Appendix A. Multi-angle case 


N 

Lo = ^ Dg • B , (15) 

N 

L = ^ Va(Da • B) , (16) 

a^l 
q -3 

where Dq, • B ~ — P^. The equations for the different 

Fourier modes [as in Eq. (11) of the main text] are 

a,n(^) ~ a.n F X Pq, 

dz 

+C50 N 

F ^ ^ ^ ^ (l ^j3,n—j ^ ■ 

3^-00 / 3^1 

(17) 


In the model considered in our work we have taken 
neutrinos to be emitted from the boundary with only two 
emission modes L and R. In this Section we generalize 
the model to a multi-angle case. As will will see the flavor 
evolution present substantial differences already without 


Numerical examples 

In this Section we show how the presence of multi¬ 
angle effects modifies the previous single-angle results 
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FIG. 5: Multi-angle flavor evolution: Two-dimensional evolu¬ 
tion of the he flavor content ^ee in the x-z plane, and its map 
on the bottom plane assuming the translational invariance 
along the ^-direction. 



FIG. 6: Multi-angle flavor evolution: Two-dimensional evo¬ 
lution of the he flavor content gee in the x-z plane, and its 
map on the bottom plane breaking the translational invari¬ 
ance along the a;-direction. 


This example is intended to capture a fundamental fea¬ 
ture of the non-isotropic neutrino emission from a super¬ 
nova core. For this reason we consider the same case of 
the work, assuming that each point along the boundary 
emits N = 200 angular modes with emission angles -da 
equally distributed in the range [0,7r]. We checked that 
the behavior of the flavor conversions in qualitatively sim¬ 
ilar further increasing the number of emission modes. In 
Fig. [5] we consider a case in which the translational in¬ 
variance along the x-direction is assumed. We realize 
that instead of the periodic flavor conversions observed 
expected in the single-angle case one now finds a flavor 
decoherence with gee{x^ z) ~ 0.5 as soon as flavor conver¬ 
sions start at z ^ 6. This case has been widely studied 
in 0 where it has been realized that in the presence of 
a constant neutrino potential, the multi-angle decoher¬ 
ence is unavoidable in both normal and inverted mass 
hierarchy. 

We show now how this result is affected by the break¬ 
ing of the translational invariance. At this regard, in 
Fig. |6] we consider the multi-angle versions of the flavor 
evolution shown for the single-angle case in Fig. 1 of the 
main text. We realize that the map of the gee(x,z) is 
rather different with respect to the corresponding one in 
the single angle case. In particular, when flavor conver¬ 
sions start they quickly lead to a flavor equilibrium with 
geeixyz) ~ 0.5 across the plane. The flavor variations 
along the x direction are smoothed with respect to the 
single-angle case being at most ~ 20%. From this figure 
one realizes that the output of the flavor evolution for a 
constant g, is the multi-angle decoherence, as also found 
in the one-dimensional models 

In Fig. [3 we show the evolution of the lepton numeber 
Lq in the x-z plane. With respect to the corresponding 
single-angle case [Fig. 2 of the main text] we realize that 
the spatial variations in Lq are ~ 0.3. Then, they are 
strongly reduced with respect to the single-angle case. 
However, also in this case we can have the formation of 
domain of lepton number with opposite sign. 

Finally, in Fig. [5] we compare the evolution of the dif¬ 
ferent Fourier modes |Pn(z)| (summed over the emission 
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FIG. 7: Multi-angle flavor evolution: Map of the lepton num¬ 
ber Lq. 
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FIG. 8: Multi-angle flavor evolution: Gontour plot of the first 
600 Fourier mode |Pn| (in logarithmic scale) in the plane n-z. 


angles) in the plane of n-z, for the case discussed above. 
We realize that the number of excited Fourier modes is 
strongly reduced (by an order of magnitude) with re¬ 
spect to the corresponding single-angle case [see Fig. 3 
in the main text]. The reason is that when the deco¬ 
herence quickly starts the length of jPo] is significantly 
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shortened. As one realizes from the non-linear term at 
second line in Eq. (113 this effect reduces the growth of 
the higher order harmonics. Moreover, the multi-angle 
kernel in the higher order terms leads itself to a suppres¬ 
sion of the higher order harmonics. This reduced growth 
of the higher order modes causes a strong suppression in 
the spatial variations of the flavor content. 


Appendix B. Analogy with non-linear flnid 
instabilities 


As discussed before, the transition of the neutrino gas 
from a coherent to incoherent behavior shows an intrigu¬ 
ing analogy with a streaming flow changing from lami¬ 
nar to turbulent regime. In this Section we discuss more 
closely the physical ground of this similarity, closely fol¬ 
lowing the presentation of [23l - l^ concerning non-linear 
instabilities of a fluid. The fluid motion is governed by 
the following equations [2^ 


V-u = 0 , 

9u ^ Vp 

-f u • Vu =- 

ot p 




(18) 


where the first equation is the continuity equation, while 
the second are the incompressible Navier-Stokes momen¬ 
tum equations, where u = (m, n, w) are the velocity com¬ 
ponents in the (x, y, z) directions respectively, p is the 
pressure, t is the time and v the viscosity. The non¬ 
linear term u • Vu is the one responsible for the turbu¬ 
lent behavior of a fluid. Eq. (fTHll admits a steady solution 
u= U{y), V = w = 0 for properly defined boundary con¬ 
ditions. In order to study the stability of this steady 
solution, one superimposes a perturbation on the basic 
flow, obtaining as disturbed velocity 


u= {u,v,'w) = {U{y) +u',v',w') . (19) 

Substituting Eq. (I19|) into Eq. m one obtains a system 
of equations for the disturbance. One can express the 
disturbance as Fourier integral over all possible continue 
wavenumbers 

r+co 

u'(x,y,t)=/ u(k,y,t)e**"'’"dk , (20) 


where x = {x,z) and k = {kx,kz). One can search for 
solutions of Eq. (fT^ expressing the Fourier components 
as 

OO 

u'(x,y,t) = ^ A„(k,<)U^(k,y) , (21) 

m—1 


where \Jm is the eigenfunction corresponding to the 
m-th eigenfrequency uim, and Am is a time-dependent 
amplitude-density function. One obtains the following 
equations for the amplitude-density functions 

dA 

+ iUJmAm -E E /(k,m,mi, 7712 , <) , ( 22 ) 

mi — 11712—2 


where 


/(k,m,mi,m2,t) = 



6(k,ki,m,mi,m2)Amj(ki,t)Am2(k- ki,t)dki . 


(23) 


The amplitude of the waves is determined by the linear 
and non-linear terms. Laminar flow contains only one 
wave, i.e. the mean flow. The linear term in Eq. (1221) 
represents the mean-flow convection of the velocity dis¬ 
turbance. The non-linear term [Eq. (1231) ] involves energy 
transfer among different waves satisfying resonance con¬ 
ditions so that the excited wavenumbers are discrete. All 
waves can resonantly interact with the mean flow. 


After this discussion about non-linear behavior of fluid 
it is apparent that the strategy used by us to follow the 
evolution of the polarization vectors P in flavor space 
is analogous to the one used to characterize the veloc¬ 
ity field V of a flow. Indeed, in both the cases the use of 
Fourier transform allows one to pass from a partial differ¬ 
ential equation problem to a tower of ordinary differential 
equations for the different Fourier modes. In both cases 
one is interested to study the evolution of disturbances 
with respect to a steady one-dimensional behavior. The 
amplitude of the excited modes is determined by linear 
and non-linear term. Neglecting the non-linear term one 
can apply a stability analysis to the linearized equations 
of motion as done in [2l| for the flavor problem and in 
for the Navier-Stokes equations. What we have done in 
the current work is to study the flavor dynamics in the 
non-linear regime, characterized by a convolution term 
among different harmonics. This term [Eq. (11) of the 
main text] has a structure similar to Eq. (1231) for the case 
of a fluid. The similar mathematical form of the prob¬ 
lem explains the intriguing analogy among the transition 
between laminar and turbulent regime of a fluid and be¬ 
tween coherent and incoherent behavior of the neutrino 
gas. Indeed as for a laminar fluid only one wave ex¬ 
ists, for the coherent flavor conversions only the n = 0 
Fourier mode evolves. Moreover, the interaction among 
the different harmonics that are rapidly excited breaks 
the coherent behavior of the flavor conversions for neu¬ 
trinos and leads to the passage into the turbulent regime 
for a fluid. 

We believe that this nice analogy would be further de¬ 
veloped to get a deeper understanding of the neutrino 
flavor evolution. Moreover, one would also benefit of the 
huge literature developed to study the non-linear Navier- 
Stokes equations. 
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